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Abstract 

The nickelate compound Y2BaNiOs is a spin-1 Haldane-gap antifer- 
romagnet. The compound is doped with holes on replacing the off-chain 
Y 3+ ions by Ca 2+ ions. Inelastic neutron scattering (INS) experiments 
reveal the existence of sub-gap states on doping. A recent INS experiment 
provides evidence for an incommensurate double-peaked structure factor 
S(q) corresponding to the sub-gap states. In this paper, we formulate a 
microscopic theory for the origin of the incommensurate peaks. 

P.A.C.S. Nos.: 75.10 Jm, 71.27.+a 



1. Introduction 

Hole-doped quantum spin systems exhibit a variety of novel phenomena as 
a function of the dopant concentration. The cuprates, which become high- 
temperature superconductors on doping, are the prime examples of such sys- 
tems. The cuprates have a rich T versus x phase diagram, where T is the 
temperature and x the dopant concentration 1-3 . In the undoped state, the 
compounds exhibit antiferromagnetic (AFM) long range order (LRO) below the 
Neel temperature Tn- The dominant electronic and magnetic properties of the 
cuprates are associated with the copper-oxide (CUO2) planes. On the introduc- 
tion of a few percent of holes, the AFM LRO is rapidly destroyed leaving behind 
a spin-disordered state. The hole-doped CuC>2 plane is an example of a quantum 
spin liquid (QSL) in 2d. Simpler AFM doped spin systems include the two-chain 
ladder compounds LaCuO^.^ and Sri^_ x Ca x Cu240Ai b and the spin-1, linear 
chain AFM compound Y2- x Ca x BaNiO^ . These quasi-ld systems exhibit phe- 
nomena some of which are similar to those observed in the cuprate systems. For 
example, the doped ladder compound Sri4_ x Ca x Cu2404i 5 ' 7 becomes a super- 
conductor under high pressure (T c ~ 12K at a pressure of 3 Gpa). As in the 
cuprate systems, holes form bound pairs in the superconducting state. 
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The doped spin-1 compound Y2- x Ca x BaNiO^ provides an example of a QSL 
in Id. The parent compound Y2BaNiO§ is a charge transfer insulator containing 
Ni 2+ (S = 1) chains. The ground state of the system is spin-disordered and the 
spin excitation spectrum is separated from the ground state by an energy gap, 
the so-called Haldane gap 8 . The compound is doped with holes on replacing the 
off-chain Y 3 + ions by Ca 2+ ions. The holes mostly appear in oxygen orbitals 
along the NiO chains. No metal-insulator transition takes place but the DC- 
resistivity pdc falls by several orders of magnitude 6 . This indicates that the 
holes are not fully mobile but delocalized over several lattice spacings. Inelastic 
neutron scattering (INS) experiments reveal the existence of new states within 
the Haldane gap. Several studies have been carried out so far to explain the 
origin of the sub-gap states 9-14 . 

A very recent neutron scattering experiment provides evidence for an incom- 
mensurate double-peaked structure factor S(q) for the sub-gap states 15 . The 
INS intensity is proportional to the structure factor. For the pure compound, 
S(q) near the gap energy of 9 mev has a single peak at the wave vector q = tt, 
indicative of AFM correlations. For the doped compound, S(q) has an incom- 
mensurate double-peaked structure factor, for energy transfer u> ~3-7 mev, with 
the peaks located at q — tt ± 5q. The shift 5q is found to have a very weak de- 
pendence on the impurity concentration x for x in the range xG[0.04 , 0.14]. 
Evidence of incommensurate peaks has also been obtained in the underdoped 
metallic cuprates 1 . The peaks are four in number and occur at (tt ± Sq, tt) and 
(tt,tt ± Sq ). The crucial difference from the nickelate compound is that Sq is 
proportional to the dopant concentration x. The incommensurability has been 
ascribed to inhomogeneous spin and charge ordering in the form of stripes. Re- 
cently, Malvezzi and Dagotto 16 have provided an explanation for the origin of 
spin incommensurability in the hole-doped S =1 nickelate compound. They have 
studied a two-orbital model using the Density Matrix Renormalization Group 
(DMRG) method. They have shown that a mobile hole generates AFM correla- 
tions between the spins located on both sides of the hole and this is responsible 
for the spin incommensurability seen in experiments. Xu et al 15 have given a 
different explanation for the origin of incommensurability. The holes doped into 
the QSL ground state of the S = 1 chain are located on the oxygen orbitals and 
carry spin. They induce an effective ferromagnetic (FM) interaction between 
the Ni spins on both sides. The incommensurate peaks arise because of the spin 
density modulations developed around the holes with the size of the droplets 
controlled by the correlation length of the QSL. In this paper, we provide a 
microscopic theory of the origin of spin incommensuration in keeping with the 
suggestions of Xu et al. 

2. Microscopic model 

Various microscopic models have been proposed to explain the presence of sub- 
gap states 9-14 . The relevant orbitals of Ni 2+ in the Ni-0 chain are 3d 3z 2_ r 2 and 
3d x 2_ y 2 13 . An electron in the latter orbital is almost localised while the 3d 3z 2_ r 2 
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orbital has a finite overlap with the 2p z orbital of the oxygen ion. In the undoped 
state, each Ni orbital is occupied by a single electron and the oxygen orbital 
is filled up. The spins on the Ni orbitals have a FM alignment, due to strong 
Hund's rule coupling, giving rise to total spin S = 1. The nearest-neighbour 
(NN) spins S = 1 interact via AFM superexchange interaction mediated through 
the intermediate oxygen ions. On doping, holes are predominantly introduced 
in the oxygen orbitals. The holes have effective spin-1/2 as each 2p z orbital 
containing a hole is occupied by a single electron (Fig. 1). 

In the undoped state, we have a chain of S = 1 spins. The ground state 
is spin-disordered and, as predicted by Haldane 8 , there is a gap in the exci- 
tation spectrum. Affleck, Kennedy, Lieb and Tasaki (AKLT) have provided 
a physical picture for the ground state. Each spin 1 can be considered to be 
a symmetric combination of two spin-l/2's. Each spin-1/2 forms a singlet or 
valence bond (VB), ^(T-l — It); with a NN spin-1/2. The ground state is 
described as a valence bond solid (VBS) in which a singlet or VB is present 
between each NN pair of spins. Actually, AKLT wrote down the Hamiltonian 
for which the VBS state is the exact ground state. The Hamiltonian is a sum 
over projection operators onto spin 2 for successive pairs of S = 1 spins and can 
be written as 



Since, in the VBS state, a singlet exists between two NN spin-l/2's, the total 
spin of two NN spin-l's can never be 2. Thus, the projection operator onto spin 
2 acting on a NN pair of spins in the VBS state gives zero. The VBS state is the 
exact ground state of Haklt with eigenvalue zero. The VBS state also provides 
the correct physical picture for the ground state of the S =1 chain interacting 
via the usual AFM Heisenberg Hamiltonian 



This has been verified in a number of experiments. 

Excitations from the VBS state are triplets (a VB is replaced by a triplet) 
which propagate along the chain 18 . The triplet excitation spectrum is sepa- 
rated from the VBS state by the Haldane gap, of the order of J. The triplet 
excitation spectrum has been observed experimentally in Y2BaNiO§. The exci- 
tation differs from a spin wave in that it propagates in a spin background with 
no conventional magnetic order. The Haldane gap A has a magnitude A ~ 9 
mev 15 . The VBS state is not, however, a simple paramagnet but has quantum 
coherence in the form of a long-range string order 19 . On doping the nickelate 
compound with holes, the triplet excitation spectrum has a dispersion relation 
similar to that in the pure Y2BaNiOs. In addition, however, new states appear 
in the Haldane gap. Near 3 to 7 mev, the magnetic scattering is found to be 
the most intense along two vertical lines, displaced symmetrically from ir. 



Haklt 



J]rfe4 + i(5l4) 2 + ! 



(1) 




(2) 



i=l 
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The oxygen hole acts as a spin-1/2 impurity and also modifies the strength 
of the superexchange interaction between neighbouring Ni spins. Sorensen and 
Affleck 10 have considered S = 1 spin chains with site/bond impurities to explain 
the origin of sub-gap states. The holes are, however, not localized but mobile 
over several lattice spacings. Various models of fully mobile S = 1/2 holes 
interacting with S = 1 spins have been proposed so far 11 ! 13 . 20 . 21 . The first such 
model is that of Penc and Shiba 11 . The effective Hailtonian is given by 



H = H + Hji + H~+ Hj (3) 



where 



jj;[i4i-^.4i) J ]. (4) 



The special cases f3 = and (3 = — -| correspond to the Heisenberg and AKLT 
Hamiltonians respectively. The Hamiltonian 

Hji = 2J'^2(T>i. ~o* l+1 /2 +6^+1.0^+1/2) (5) 

i 

describes the effect of S = 1/2 impurities (holes). An impurity at site i+1/2 
changes the J coupling to J\ between the S =1 spins on both sides of the impu- 
rity in the Hamiltonian H (Eq.(4)). Also, there is an interaction between the 
impurity spin a and the adjacent S = 1 spins, -ff-and H ~ are the Hamiltonians 
describing the effective hopping of the holes on the sites. 

H 7 =tJ2Pi (6) 

i 

and 

Hj= 2J^P i (^ i .lf i + 1/2 +l? i _ 1/2 .'3 i ). (7) 

i 

The operator Pi exchanges the occupations of the sites i +1/2 and i— 1/2 . If 
the hole is located on site i+1/2, Pi shifts it to i-1/2 provided that site is empty. 
Two S = 1/2 objects cannot occupy the same site as it costs a lot of energy. 
The same type of effective Hamiltonian has been used earlier to describe doped 
cuprate systems in a two-band (copper and oxygen orbitals) scenario 22-24 . In 
the cuprates too, the added holes go to oxygen sites and each hole has an effective 
spin-1/2 which interacts with the spin-l/2's of the Cu 2+ ions on neighbouring 
sites. The major differences between the nickelate and cuprate compounds are: 
Ni 2+ ion carries spin 1 as opposed to the spin-1/2 of the Cu 2+ ion and in the 
undoped state the latter compound has AFM LRO. 

For the doped cuprate systems, Aharony et al 25 have argued that the holes 
generate an effective FM exchange interaction between the Cu ions on neigh- 
bouring sites. Consider a pair of Cu ions and an intermediate hole. The in- 
teraction is described by a Hamiltonian of the type given in Eq.(5) except that 
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the spins Si have a magnitude 1/2. The exchange interaction strength J' is 
expected to be larger than the interaction strength J. The spins Si and S^+i 
prefer to be parallel irrespective of the sign of J'. The spin-1/2 of the hole 
along with the two spin-l/2's of the adjacent Cu-ions constitute a spin bubble 
or polaron. For the doped Y2BaNiO§ compound, a similar scenario holds true. 
The hole with a spin-1/2 provides a break in the VBS state of the S =1 chain. 
The two Ni-ions which are NNs of the hole have unpaired spin-l/2's. In the 
absence of the hole, these two spin-l/2's form a singlet in the VBS state. In the 
presence of the hole, a spin bubble consisting of three spin-l/2's (the central 
spin belonging to the oxygen hole) is obtained. The bubbles are embedded in 
the original VBS state. 



3. Structure factor S(q) 

The differential scattering cross-section in an INS experiment is proportional to 
the dynamical structure factor S aa (q,u>), (a = x,y or z) given by 26 



S aa (q,cj) = Yl Y. eX P^ - r i)> < 5M I S ? \ S ' M ') < 5 ' M 'I S J \ SM ) 

S',M' i,j 

5{lj + E{S,M)- E(S',M')). (8) 

We assume a Id spin system at T = 0. The magnetic ions are located at the sites 
i and j with the associated magnetic form factors taken to be unity. The spin 
state \SM) is the ground state where S is the total spin of the state and M the z 
component of the total spin. The ground state energy is given by E(S,M). The 
state \S'M') is an eigenstate of the system with energy E(S' , M'). The energy 
and momentum transfers in the scattering process are u> and q respectively. 

For the doped Y2BaNiO$ compound, the spin states are the VBS states 
with spin bubbles nucleated around the doped holes. Let us consider the case of 
a single static hole. The total spin of the state is equal to the total spin of the 
bubble since the rest of the VBS state is a singlet. Since the bubble consists of 
three spin-l/2's, the possible spin states correspond to total spin S = 3/2 and 
1/2 (two representations). The possible spin states \SM) are (only the bubble 
spin configurations are shown): 



3 3 
2' 2 
3 1 
2' 2 
3 1 
2' _ 2 
3 3 
2' _ 2 



ITTT) 



= _|tiT + m + TTT) 
= -j= ITTT + TTT + TTT) 
= ITTT) 



(9) 
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1 1 

2' 2 



^=|2UT-m-m> 



i 

2' 



V6 



|2m-TU-UT> 



(10) 



^) = ^im-m) 



i'-i) = ^i IUT - TU> 



(ii) 



The '"'-sign in Eq.(ll) indicates a different S = 1/2 representation. In the 
states shown in Eqs.(9)-(10), the first and the third spins (the Ni spin-l/2's) 
are in a S = 1 spin configuration. This spin 1 adds to the spin 1/2 of the hole 
giving rise to a resultant spin 3/2 or 1/2. In all these states, there is thus an 
effective FM interaction between the Ni spins. In the states shown in Eq. (11), 
the Ni spin-l/2's form a singlet and the hole spin is free. 

Penc and Shiba 11 have calculated the hole excitation spectrum assuming a 
fully mobile hole. The trial wave function of the L-site periodic chain with one 
hole is given by 



\SMk) = ^=J2 eikj 

j 



SM,j + - 



(12) 



where 



\SM,j + 5) denotes the state with the hole located at site j + |. With 
Hq in Eq.(4) given by the AKLT Hamiltonian, the lowest hole eigenstate of H 
Eq.(3)) is given by an appropriate linear combination of the states j^^fe) and 

\\k^- An excited state is obtained from a different linear combination of the 

states. The full hole spectrum, in the subspace of states described by Eqs. (9)- 
(11), consists of one fourfold and two twofold degenerate bands. The dynamical 
structure factor has non-zero weight inside the Haldane gap due to scatter- 
ing between the hole states. Penc and Shiba's microscopic model Hamiltonian 
has been rederived by Batista et al 14 starting from a multiband Hamiltonian 
containing the relevant Ni and O orbitals. They have further estimated the 
magnitudes of the various parameters in the Hamiltonian on the basis of small 
cluster calculations. 

We return to our consideration of a static hole. From Eq.(8), the selection 
rules for INS transitions are 



AS = 0, ±1, AM = 0,±1. 



(13) 



We consider a = z so that AM = . The q-dependence of the INS intensity is 
given by 
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is'M'(q) = 



(S'M'^e^'Sf^M) 



(14) 



Is'M'(q) corresponds to a scattering transition from the state \SM) to the state 
\S'M') with an energy transfer of uj. The spin expectation values in Eq.(14) 
can be calculated in a straightforward manner using the method of matrix 
products 2 . In the undoped state, the Hamiltonian H in Eq.(3) reduces 
to Hq . We consider Hq to be the AKLT Hamiltonian for which the VBS state 
is the exact ground state. In the matrix product representation, the VBS state 
is given by 

\tpVBs) = Tr{g_ L <g> .... <g> g Q <g> .... <g> g L ). (15) 

The Id lattice consists of 2L + 1 sites and gi denotes a 2 x 2 matrix at the ith 
site with single spin states as elements. 



W -|0) -V2|l) 
9 V3^V2|1) |0) 



(16) 



where |1) , |0) , |l) are the spin-1 states with S z = +1,0 and —1 respectively. 
Expectation values in the matrix product states can be calculated most con- 
veniently by the transfer-matrix method. We describe the major steps of the 
calculation in the following. Details of the theory of matrix product represen- 
tations can be obtained from Refs. [27-30]. 

The norm of a wave function in the matrix product representation, e.g., 
(ipVBS | i>VBs) can be written as 



(ipVBS I ipvBs) = ^2 

(17) 

{n a ,m a } 



We define a 4 x 4 transfer matrix G at any lattice site by 



G„ 



G 



(n 1 m 1 ),(n 2 m 2 ) 9n 1 n 2 9m 1 m 2 



The ordering of multi-indices is given by 

H = 1,2,3,4 <-> (11), (12), (21), (22). 
With the choice of g given in Eq.(16), the transfer matrix G is 



G 



/ 1 2 \ 
1 0-100 

3 0-10 

\ 2 1 / 

The eigenvalues and eigenvectors of G are 



(18) 
(19) 

(20) 
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Ai = 1, A2 = A3 = A4 = — — 



(21) 
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The normalisation (17) is determined by the largest eigenvalue for L — > 00 (the 
size of the chain is 2L+1): 

(ipvss I tPvbs) - Af +1 = 1. (23) 

Consider the ground state matrix element of a spin operator P at site r. One 
uses the notation 

Z(P)^ 2 Z(P) 

(n 1 m 1 ),(n 2 m 2 ) ~ 9n 1 n 2 ^ > 9mim 2 (24) 

with Z(1)=G. We obtain 



(P) = (^vbs I ^VBsT 1 0>vbs I P I 4>vbs) = {TrG 2L+1 )- 1 TrZ(P)G 2L . 

(25) 

On taking the limit L — > oo , we obtain 

(P) = \^(e 1 \Z(P)\e 1 ). (26) 

For 2-site correlations of operators A\ and B r at sites 1 and r respectively, one 
obtains 

(Aifl r ) = Tr{G 2L+1 )- 1 TrZ(A)G r - 2 Z{B)G 2L+1 - r . (27) 
For L — > oo, the expectation value reduces to 

(^fl,.) = 2 A^ 2 (^Y (ei| \e n ) (e n \ Z{B) | ei > (28) 

where A„ and |e„) are the eigenvalues and eigenvectors of the transfer matrix G. 
In the same way, one can calculate off-diagonal expectation values of the type 
(S'M'\Ai\SM). 

The VBS state with a single embedded hole has different matrix product 
representations depending on the configuration of the three-spin bubble centred 
around the hole. The possible bubble spin configurations are given in Eqs. 
(9)-(ll). For example, the states |||) and |||) are of the form 

Tr{g_ L ® ®g-i®g h ®g ®....®g L ) (29) 
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with the hole located in between the sites -1 and 0. The corresponding matrix 
g h is given by 



for the HI) state and 



9h 





V2|T> o 



-IT) o 
-2 II) IT) 



(30) 



(31) 



for the 



ill' 



.-. 2 / StcltG. 

We have calculated Is'M'{q) in Eq.(14) by taking the ground state \SM) to 
be the \ \\) state and the state \S'M') to be a state which satifies the selection 
rule AS — 0, ±1 and AM = . For example, the state \S'M') can be any one 

of the states |||), \\\) and \r^J- In the first two cases, Is'M'(q) is calculated 

as (apart from numerical prefactors ) 



is'M'(q) = 



(i 



.K 



)cos(!) 



(32) 



coshK + cosq 

where K~ x = j s the spin correlation length in the VBS state. For the state 

H2) ' 



9h 



IT) 

-II) 



o 

-IT) 



(33) 



Also, i , j in Eq.(8) denote the sites at which the Ni spins are located, i.e., the 
hole spin is not explicitly taken into account. When \S'M') is 



9h 



IT) o 
o |T) 



and Is'M'(q) is given by (apart from a numerical prefactor) 

is'M>(q) = 



(i 



sin(i) 



(34) 



(35) 



coshK + cosq 

Xu et al 15 have suggested a form (Eq.(32)) for the structure factor of the ground 
state with a single FM bond inserted within an infinite chain. The FM bond 
occurs between the sites -1 and 0. In the matrix product representation, such a 
state is given by 



where 



\iPfm) = Tr{g_ L ® ® g a _ x <g> g 



a _ ( V2|l) 



>9l) 



-|0) 



(36) 



(37) 
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Calculation of S(q) with \S'M') = \SM) = \ipFNi) reproduces the form in 
Eq.(32). This, however, corresponds to elastic neutron scattering whereas the 
incommensurate signal arises from inelastic neutron scattering. The intensity 
Is'M'(l) (Eq.(32)) has two peaks symmetrically displaced from ir. The widths 
and incommensurability are of the order of K. A good fit to the experimental 
data has been obtained by treating K~ x as a parameter. For dopant concen- 
trations x = 0.04, 0.095 and 0.14, the best fits to the experimental data are 
obtained for R- 1 = 8.1 ± 0.2, 7.3 ± 0.2 and 7.2 ± 0.5 respectively. These values 
are considerably higher than the VBS value of K^ 1 = and are closer to 
the estimate of the correlation length for the Heisenberg Hamiltonian. The use 
of the AKLT Hamiltonian has, however, enabled us to formulate a microscopic 
theory of the structure factor, based on the matrix product representation. The 
double-peaked incommensurate structure is correctly reproduced in the theory. 
The intensity Is'M 1 (q) , however, has nodes at q = (2n+ l)ir which are not seen 
experimentally. Xu et al 15 have suggested that finite hole densities as well as an 
explicit consideration of the hole spin distributed over 1 lattice sites centred on 
the FM bond, gives rise to a non-zero intensity at q = (2n + 1)tt. Microscopic 
calculations along these lines are, however, yet to be carried out. 

The structure factor S(q) of INS is obtained on integrating the dynamical 
structure factor S(q, lo) in Eq.(8) over energy ui and is given by the sum over the 
intensities Is'M'iqYs for the three possible states \S'M'). Two of the intensities 
is proportional to cos 2 | while the third is proportional to sin |. The first two 
intensities have a minimum at q = ir and the last intensity has a maximum 
at q = 7T. The resultant intensity for K = ln3 <~ 1.1 has a double-peaked 
incommensurate structure but the intensity at q = ir has a magnitude very close 
to that of the peak intensities. For K = 1.2, (correlation length less than NN 
Ni-Ni distance but more than the distance separating NN Ni and oxygen ions), 
the structure factor S(q) vs. q is shown in Fig. 2. The double-peaked structure 
is, however, very sensitive to the value of K and disappears for K < 1.0. For 
exchange interactions of the Heisenberg-type, K is around 0.17. With a ground 
state \S M) = \ \\) +x (x < 1), the structure factor S(q) contains terms 

proportional to cos 2 | as well as Ssin 2 ^ (S < 1). With the appropriate choice 
of S, a double-peaked incommensurate structure factor S(q) is obtained for a 
wide range of values of K. Fig. 3 shows S(q) vs. q for S — 0.35 and K = ln3. 

Eigenfunctions which are linear combinations of the states \ \\) and |^ are 
obtained when mobile holes are considered. In fact, as mentioned earlier, for 
fully mobile holes , the lowest two hole bands correspond to linear combinations 



of the states and \\kj with \SMk) defined as in Eq.(18). Holes in the 

doped nickelate compound are not, however, fully mobile. Calculations, taking 
into account the limited mobility of a hole as well as an explicit consideration of 
the hole spin in calculating S(q), are in progress and the results will be reported 
elsewhere. 
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4. Conclusion 



In this paper, we have provided a microscopic theory for the double-peaked 
structure of the INS intensity for the doped S = 1 AFM compound YzBaNiO^. 
The spin-l's are assumed to interact via the AKLT Hamiltonian for which the 
VBS state is known to be the exact ground state. The introduction of holes 
in the system nucleates spin bubbles around the holes due to the strong AFM 
exchange interactions between the hole spins and NN free spin-l/2's of the VBS 
state. The scattering between the hole states gives rise to states within the 
Haldane gap with an incommensurate structure factor. The incommensurate 
peaks occurring in the underdoped metallic cuprate systems have been ascribed 
to the formation of spin and charge inhomogeneities in the form of stripes. As 
a consequence of this, the shift Sq of the location of an incommensurate peak 
from 7r is proportional to the dopant concentration x. For the doped nickelate 
system, 5q has a very weak dependence on x for x in the range 0.04-0.14 and 
no experimental evidence of stripes has been obtained. Incommensuration in 
the doped cuprate systems occurs in the metallic state when holes are fully 
mobile. In the doped cuprates, the spin-spin correlations have a longer range 
than in the case of the doped nickelate system. This gives rise to greater effective 
interactions between the spin bubbles in the former case. 

Malvezzi and Dagotto 16 have considered a two-orbital model for the doped 
Y2BaNiO$ compound and have shown that spin incommensurability is obtained 
in a wide range of densities and couplings. AFM correlations are dynamically 
generated across the holes to facilitate hole movement. The same type of corre- 
lations occur in the Id Hubbard model and the 2d extended t-J model. In the 
two-orbital model, only the two Ni orbitals have been considered and oxygen 
orbitals have not been explicitly included in the Hamiltonian. The spin of the 
hole on the oxygen ion forms a singlet with the spin- 1/2 in the adjacent Ni 
orbital (each of the two Ni orbitals contains an electron with spin-1/2) giving 
rise to an effective Zhang-Rice (ZR) doublet. This is analogous to a ZR singlet 
in the case of doped cuprate systems 2 . More detailed experiments are required 
to distinguish between the two mechanisms, suggested so far, on the origin of 
spin incommensuration in Y2- x Ca x BaNiO§. 

Several earlier studies on doped cuprate systems have explored the possibility 
of the binding of two spin bubbles in the antiferromagnetically interacting spin 
background. Zhang and Arovas 31 have considered the propagation of S = 
holes in a VBS state and have found evidence for the binding of holes. The 
propagating S = hole corresponds to a simultaneous hopping of a pair of 
electrons. The possibility for the binding of a pair of spin bubbles, centred 
around holes, in a VBS-type spin background has not been explored as yet. 
Experimental evidence of the binding of holes in the doped nickelate compound 
does not exist. On the other hand, evidence for the binding of holes has been 
obtained in the case of a doped ladder compound 3,5 . Several numerical studies 
provide evidence for the binding of holes in a two-chain t-J ladder model 3 ' 32 . 
The binding of two holes has been demonstrated exactly and analytically in a 
specially constructed t-J ladder model 33 . In the undoped case, this model, in 
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an appropriate parameter regime, has a ground state identical to that of a S 
= 1 AFM Heisenberg spin chain 34 . The exact results on the binding of two 
holes hold true in a different parameter regime. The connection between doped 
spin-1 and ladder systems is not clearly established as yet. In summary, the 
experimental results on the S = 1 doped nickelate compound provide examples 
of the rich phenomena observed in other spin systems. Comparative studies of 
such systems are needed to highlight the common origins, if any, of some of 
these phenomena. 
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Figure Captions 

Fig. 1. Schematic diagram of the Ni-0 chain with the second being occupied 
by a single hole. 

Fig. 2. Structure factor S(q) vs. q for K = 1.2 

Fig. 3. Structure factor S(q) vs. q for K = ln3 and S = 0.35 
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